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Abstract. We investigate wave amplification through the 
Electron Cyclotron Maser mechanism. We derive a semi- 
analytic approximation formula giving the frequencies at 
which the absorption coefficient is negative. The coeffi- 
cients still need to be computed to obtain the largest, and 
therefore the dominant, coefficient. 
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1. Introduction 

The observation of millisecond microwave spikes from the 
Sun has been interpreted as gyrosynchrotron maser emis- 
sion since the late 1970's ( p"l] ; ||). The spikes exhibit a dura- 
tion of a few milliseconds to some tens of milliseconds, and 
a narrow bandwidth (w 1% or a few MHz). The charac- 
teristic high brightness temperature (10 15-18 K) deduced 
for the spikes strongly suggests maser action (|||). 

Several papers have explored the possibilities of cy- 
clotron maser. Melrose and Dulk (^) have used the semi- 
relativistic approximation to derive approximate formulae 
for the frequency of the largest growth rate, and for the 
growth rate. Winglee ( |l0| ) approximated the effects of the 
ambient plasma temperature on the maser emission. As- 
chwanden (|l|, |^) followed the diffusion of electrons into the 
loss-cone as a result of the maser emission, and therefore 
the self-consistent closure of the loss-cone as the energy 
is transferred to the maser. Kuncic and Robinson @ per- 
formed ray-tracing in a model loop with dipole field, and 
concluded that maser emission can escape from lower lev- 
els. 

Some papers have also attempted to derive analytical 
approximate formulae for the frequencies of the maser. 
Dulk and Melrose (||) concluded that the maser emission 
is at, or near, 90 degrees to the magnetic field and the 
maximum growth rate they find is very near the cyclotron 
frequency, or at the second harmonic. Hewitt and Melrose 



(^) derive the conditions for solutions and the number of 
solutions, but not the frequencies themselves. 

All previous studies concentrated on the growth rate 
of the amplified wave, we use the absorption coefficient 
approach of Ramaty (Q) and calculate the frequencies at 
which the absorption coefficient is negative. 

2. Gyrosynchrotron Absorption 

In studies of microwave emission from solar flares it is 
customary to assume that the plasma is composed a non- 
thermal population of energetic elecrons, and a thermal 
population of much higher density, but much lower energy. 
The thermal plasma is assumed to determines the propa- 
gation of waves. This plasma can be described by the cold 
plasma approximation using magnetoionic theory (]?]). 

ECM studies usually compute the growth rate, while 
we follow Ramaty (^|) in using the absorption coefficient, 
but for the purposes of this article the detailed equations 
are unimportant. Both the growth rate equation and the 
absorption coefficient equation include a 5-function from 
which the important resonance condition is derived 
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where lub = eB/m e /c is the cyclotron frequency, to is the 
emission frequency, s is an integer, 9 is the angle of the 
wave vector k to the magnetic field, <p is the angle between 
the velocity of the electron and the magnetic field (the 
pitch angle), and 7 and (3 are the usual Lorentz factor 
and velocity. The important factor in equation [l] is n± 
the refraction index for the Ordinary Mode (OM or '+'), 
and the eXtra-Ordinary mode (XO or '-'). The refraction 
index is a function of the ambient density, of the cyclotron 
frequency, of the emission frequency, and of the angle to 
the magnetic field 
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Where u p is the plasma frequency and the rest are as in 
equation |TJ. 

If 9 7^ ir/2 the 5 function is transformed to 



d I v — — : — — n±v(3cos((j)) cos(9) ) = 
V 7 
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This transformation gives us the pitch angle cf> s of an elec- 
tron emitting in direction 9 a wave of frequency v, where 
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The absorption coefficient after the integration over 
the pitch angle is (||) 
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Where J s are Bessel functions, with the argument x s = 
{v/vB)n±^j3sm{6) sin(<^)). The cold plasma parameter T± 
is the ratio of the polarization ellipse axis, which describes 
the transverse part of the polarization vector, and the pa- 
rameter L± describes the longitudinal part of the polar- 
ization vector. Equations for T± and L± can be found in 
Melrose ([?]). The range of the summation over the Bessel 
functions is 



S12 = — (l±ra±/?cos(0)) 
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And u(-y)g((f>) is the distribution function of the emitting 
electrons. 

We assume the distribution function in equation ^| can 
be put in the form 

f(p)d 3 P = 27ra(7)g(0)d7d(cos(<£)) 

And take for the distribution function of the non-thermal 
electrons a power-law 



u ( 7 ) =A-(j-iy 



(7) 



which is normalized to n^ ot the number density of the en- 
ergetic electrons. For the pitch angle dependence in the 
energetic distribution we use an ideal loss-cone which de- 
scribes an isotropic distribution for pitch angles above 



some a, no electrons below another angle a — Aa, and 
a linear in cos(4>) decrease within the Aa interval (ra), or 



g(<t>) 



A, cos(0) < cos(a) 

Ax Zt-tt-coslt] , cos(a) < cos(0) < (g) 

cos{a — Aa) 
0, cos(a — Aa) < cos(<p) 



Where A is a constant which depends on the loss cone 
parameters, and for an idealized loss cone as above is A = 
1/tv/(2 + cos(a) + cos(a — Aa)). With this constant the 
integral is normalized so that 2tt f , g((j))d(cos <j>) = 1. 

3. Conditions for Maser 

The Electron-Cyclotron-Maser occurs when the absorp- 
tion coefficient is negative. The resulting amplification of 
the radiation is referred to as maser. 

Inspection of equation || shows that the absorption co- 
efficient for a distribution which is isotropic in pitch angle 
can only be negative if the slope of u(7)/7 2 //3 is positive. 
This is the equivalent of an increase of the fast particles 
with energy, a distribution which is not likely to occur. 
The other possibility for amplification is when the dis- 
tribution is not isotropic in pitch angle. In that case if 
the derivative of g(4>) is large enough, and there are solu- 
tions such that the pitch angle dependent term in equa- 
tion |^ is negative, then the absorption will be negative. 
One such anisotropic distribution is the loss-cone distri- 
bution, which we use. The loss-cone distribution produces 
negative absorption coefficients and maser emission. 

4. Semi-Analytical Formulae 

The resonance condition |l| defines an ellipse on the elec- 
tron velocity plane (v\\,v±) for a given s,9, and frequency 
of emission (0). The integral over the electron distribu- 
tion function, which is necessary to calculate the absorp- 
tion coefficient, is therefore an integral along the ellipses 
of the various s harmonics. In the weakly relativistic case 
1/7 = 1 — /3 2 /2 this ellipse becomes a circle (||) - the res- 
onance circle - whose center is on the Wn-axis (the x-axis) 
at v\\ = v c , and whose radius is Vr with 



v c n± cos(#)z^ 

C SVb 




2{v - su B ) 



svb 



(9) 



(10) 



The loss-cone described by equation || divides the 
(v||,i>jj plane into three regions. Below the line of angle 
a — Aa to the uii-axis there are no electrons, above the 
line of angle a to the un-axis the distribution is isotropic 
in angle, and between the two lines the number density of 
electrons decreases with angle. 
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The frequency and angle for which the largest (in ab- 
solute value) negative absorption is expected are those for 
which the resonance circle is tangent to the outer edge of 
the loss-cone at a (||), as shown in figure [l]. Integration 
along this special circle passes only through the two lower 
regions where there are either no electrons, or the distri- 
bution changes with angle. Where there are no electrons 
the absorption is zero, and where the distribution changes 
with angle is where the integrand is most likely to be neg- 
ative, and therefore integration around the special circle 
is the most likely to give negative absorption. The special 
circle also has the longest possible path through this re- 
gion, and is expected to produce the largest (in absolute 
value) negative absorption coefficient. 

The definition of the special circle results in an equa- 
tion for a frequency which is the frequency where the 
largest absorption is expected. This frequency is (H) 
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Since even a short path through the positively contribut- 
ing region of the (vt \,v±) plane increases the positive ab- 
sorption considerably, the expectation is that there are 
negative absorption coefficients only for this special fre- 
quency, or at frequencies very close to it. 

For the fully relativistic case it is more convenient to 
consider the plane (7, cos(4>)). On this plane the loss cone 
regions are the region below the line parallel to the 7-axis 
(the IE-axis) at cos(</>) = cos(a — Aa), where the num- 
ber density is independent of the pitch angle, the region 
above the line cos(0) = cos(a), where there are no elec- 
trons, and the region between the lines, where the number 
density decreases. On the (j , cos((f>)) plane the resonance 
condition defines a curve which is given by equation |j| 

There are two possible forms to the curve given by 
equation Q as shown in figure I The first for m is for 
svb > v. For this case, as 7 approaches 1 from above, and 
assuming cos(6>) > 0, the cos(0 s ) solution goes to minus 
infinity. The second form of the curve is for svb < v, 
and for this form the cos(</> s ) solution goes to positive 
infinity as the kinetic energy approaches zero. Both curves 
go asymptotically to [n±cos(0)]~ 1 as 7 goes to infinity. 
The second form curve (for v > sv B ) undergoes a change 
in the sign of the slope at some energy. For that curve the 
value of cos(0 s ) ceases to decrease as 7 increases (negative 
slope), and starts to increase (positive slope), as figure || 
illustrates. 

For the fully relativistic case the equivalent curve to 
the special resonance circle is a curve with the second 
form, which on the (7,cos(</>)) plane is tangent to the line 
cos(0) = cos(a). This curve goes only through the two 
upper regions of the plane, where there are either no elec- 
trons, or the density is pitch angle dependent - where the 
integrand in equation || is most likely to be negative. This 
special curve also has the longest possible path through 



the negative contribution region, and therefore the absorp- 
tion coefficient derived from integration along this curve 
should be the largest negative absorption coefficient pos- 
sible. The turnover of the special curve is at energy 



7 = v/{sv B ) 



(12) 



and the frequency which gives this curve is given by the 
equation 



(13) 
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The equation [l3| is our new approximation. Our ap- 
proximation becomes identical with the Melrose and Dulk 
approximation 11 for the case cos 2 (9) cos 2 (a) << 1. 



It is, however, more general and is a good approxima- 
tion for the frequencies of negative absorption, even when 
n± cos 2 (9) cos 2 (a) is large, as the examples in section || 
demonstrate. 

Even though we utilize the approximation |l3] it is im- 
portant to note that it is also possible for a curve of 
the first form, with sug > v, to give a negative ab- 
sorption coefficient. This may occur if the curve reaches 
cos(0) = cos(a) for an energy which is close to the low- 
energy cut-off of the fast electron distribution. For this 
case most of the curve at cos(0) < cos(a) is at energies 
where there are no fast electrons. This part of the curve 
therefore has no contribution to the absorption coefficient, 
and since this is the region where the positive contribution 
comes from, the resulting absorption is negative. We used 
in our calculations a low-energy cut-off of 7 ml „ = 1.02, 
which is a common assumption. It turns out that it is 
very unlikely that equation || gives cos(0 s ) = cos(a) for 
an energy lower or close to this cut-off. The equation [l^ 
therefore gives most of the solutions, and can be used to 
find the frequencies of negative absorption. 

Equation |l3| is still an approximation, and in deriving 
it we made simplifying assumptions . For example frequen- 
cies close to the special frequency given by equation 13 also 
have negative absorption, because as long as the curve for 
some frequency has a short path through the isotropic 
region the absorption can remain negative. Equation [l^ 
also assumes a single harmonic s while the absorption 
coefficient equation [5] includes a sum over all harmonics 
from s\ to S2- The positive contribution of other harmon- 
ics moves the negative absorption towards zero, and can 
not be easily approximated. Therefore the frequency with 
largest negative absorption is usually not exactly the so- 
lution of equation O, but a frequency close to it. 



A close look at equation 13 shows that it must be 
solved numerically, because the refraction index n± is a 
function of the frequency v. The refraction index is where 
the ratio between the plasma frequency v v and the cy- 
clotron frequency vb comes into play, since it is very sen- 
sitive to this ratio for low multiples of vb ■ Our approxima- 
tion is therefore only semi-analytical, since some numerical 
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computation must be made. The computation, however, 
is very simple and can be performed very quickly. 

There are some general properties of the solutions it is 
possible to derive even without calculation of the refrac- 
tion index, and we proceed to do so in the next section. 

Finally, equation [l3| was derived with the assumption 
that cos(#) > 0, but if cos(0) < then cos(a) < as well, 
and the solution is the same. Equation |l3| can therefore 
be used for all cases. 



5. General Properties of the Solution 

The first general property of the solution given in equation 
|l3| is the requirement that v > svb ■ This may not always 
be the case (as shown in the previous section), but it turns 
out that this condition holds for the dominant (defined as 
the most strongly amplified) modes. In practical terms 
this condition sets the appearance of negative absorption 
at frequencies between svb and (s + 0.5)vb- 

Other general properties can be demonstrated by a 
graphical solution of equation |l3|. On a plane where the 
x-axis are the frequency we plot the line y = v, and the 
function 



/(") = 



SV B 



yjl — n± cos 2 (6) cos 2 (a) 



(14) 



Any frequency for which the two curves intersect is a so- 
lution of equation |l3|. 

There are two special cases for the OM and the XO 
mode where we can compare f{v) and y = v without 
explicitly calculating the refraction index. The first case 
is at the special points of the cut-off frequencies of the 
modes, where n± = 0. The second special case is for large 
frequencies, where the refraction indices are n± = 1. The 
general shape of the curve |l4| can be inferred from these 
two points. 

The OM cut-off is at v = v p , and therefore for the 
Ordinary Mode f(v p ) = svb, and if the plasma frequency 
is equal to the cyclotron frequency, or some harmonic of 
the cyclotron frequency, then v v is a solution. The XO 
mode cut-off is at the frequency 



(15) 



This cut-off is always at higher frequency than vb, and 
therefore j{v x ) < v x with s = 1. However, if v p = \[2vb 
then the cut-off is at v x = 2vb, and the second harmonic 
s = 2 is a solution. For any harmonic s, the cut-off fre- 
quency is equal to svb for a ratio (v p /vb) 2 = s 2 — s. 

The cut-off frequencies of the modes are not important 
themselves, because a wave with a zero refraction index 
can not propagate. However, the solutions can be used to 
deduce the behaviour of the curve [l4|. The line y = v and 
the curve [l4| can have one of two relations at the cut-off 
frequencies : 



— If the plasma frequency is smaller than the cyclotron 
frequency then f{v p ) = vb > v p and the curve |l4| for 
the OM starts above the y = v line (figure ^). More 
generally, for all harmonics s such that svb > v p the 
curve for the OM starts above the y = v line. 

— Conversely, if the plasma frequency is larger than some 
harmonic m of the cyclotron frequency, then the OM 
curve |l4| starts below the y = v line for s up to and 
including the harmonic m (figure ^). 

— If (v p /v B ) 2 < s 2 - s, then f(v x ) = sv B > v x , and 
the curve [l4| for the XO starts above the line y = v. 
because then v x < svb- 

— If (v p /vb) 2 > s 2 — s, then the XO curve starts below 
the line y = v. 

The behaviour of the curve [li] for frequencies larger 
than the cut-off frequencies is deduced from the refraction 
index. Both the OM and the XO mode have refraction 
indices smaller than 1 which approach 1 asymptotically as 
the frequency increases. A curve like the curve [l4| therefore 
has the shape illustrated in figure || if it starts above y = v. 
From figure || it is seen that for a starting point above the 
y = v line there is only one intersection. If the curve starts 
below the y — v line there are either two intersections, or 
no intersection at all, as is illustrated in figure |I| The 
conclusion is that for the OM there is one solution for 
negative absorption coefficients if svb > v v , and for the 
XO there is one solution if \/ s 2 — s vb > v p . There are 
either two solutions or no solution for the OM if svb < v p , 
and for the XO there are two solutions or no solution if 
\/ s 2 — s vb < v p . In practice for the cases where the curve 
|l4| starts below the y = v line, there is usually no solution, 
and only the case illustrated in figure^ is important. 

There is an additional limit we can derive from equa- 
tion |l3|. The mathematical formula derived in equation [l4| 
can be plotted as a function of frequency for any frequency. 
However, beyond some frequency the mathematical result 
f(v) is no longer relevant to finding the negative absorp- 
tion coefficient. The most important limit is derived from 
equation [TJ itself, and it is the highest frequency possi- 
ble for a given harmonic s with this equation. Clearly, 
plotting f(v) beyond this frequency is meaningless. The 
highest possible frequency is derived for the OM and the 
XO by assuming that n± = 1, and is therefore : 



sv B 



sjl — cos 2 (0) cos 2 (a) 



(16) 



For example with cos(a) = 0.8, s = 1, and emission di- 
rection cos(#) = 0.3 the upper limit on the frequencies 
usable for this set of parameters is v < 1.03vb- If we take 
cos(#) = 0.7 the limit is v < 1.207i^b. This result gives 
the range of frequencies where negative absorption is ex- 
pected. For example, cos(#) = 0.3 can only have negative 
absorption between v = vb and v = I.OZvb- 

The condition [l6] can also be used to limit the harmon- 
ics relevant for any frequency and angle of emission. For 
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example the first harmonic s = 1 is not useful for find- 
ing negative absorption coefficients at u > 2ub for any 
loss-cone with cos(a) < 0.866. It is thus possible to im- 
mediately disqualify certain harmonics from contributing 
in certain frequency ranges, and shorten the computation 
time further. 

In summary we conclude that for the OM and the XO 
mode, we expect negative absorption coefficient to appear 
for frequencies obeying the following two conditions : 

— For the OM the frequency should be derived from a 
harmonic s which has s > u p /ub- 

— For the XO mode the frequency should be derived from 
a harmonic s which has s 2 — s > {u p /ub) 2 ■ 

— For both modes the frequency must be smaller than 
v-max of equation |l6| 

The harmonic number s of a frequency can be derived 
from the condition [l(] and should be the smallest integer 
which fits the inequality 

s > —y/l - cos 2 (0) cos 2 (a) (17) 
vb 

The largest harmonic s which contributes to a frequency 
is given by equation ^, but it is reasonable to assume that 
the smallest harmonic has the largest contribution, and 
therefore we can consider it alone for the purposes of the 
approximation. 

So far we have not considered the important Z-mode 
(ZM). The Z-mode is the lower branch of the XO mode, 
which reappears at frequencies smaller than u z 

= 4 + »i + v^i + ^ 2 ) 2 - cos2 w ( is ) 

and its cut-off frequency is u x — vb- 

Emission in the Z-mode can not emerge from the 
plasma, and therefore can not be observed. Since our pur- 
pose is to derive estimates for the observable frequencies, 
approximate solutions for the Z-mode are not of great in- 
terest to us. However, the Z-mode may be the dominant 
mode and quench the maser before the other modes are 
amplified. It is therefore of interest to determine where 
possible negative frequencies can appear. 

The upper Z-mode frequency u z is bound from above 
by the upper hybrid frequency 

vuh = + (19) 

From the equation it is clear that u z decreases with in- 
creasing cos(0), and for cos(0) = 1 the limiting frequency 
is either v z = ub if u p < vb or v z = v v if u p > Vb- The 
lower frequency boundary of the ZM is always smaller than 

Up. 

The ZM refraction index is n re f > 1 for frequencies 
u > Up, is n re f = 1 at the plasma frequency, and is smaller 
than 1 for frequencies smaller than u p . 

Using the above characteristics we can determine the 
general behaviour of the curve |lj for the Z-mode, which 



is illustrated in figure [5| The curve goes to infinity for 
frequencies near u z , and goes to sub as the frequency ap- 
proaches u x — ub- 

As for the OM and the XO mode it is possible to de- 
duce the general behaviour of the solutions for the Z-mode 
from the properties of the refraction index, however, be- 
cause of the additional dependence on cos(0) the ZM prop- 
erties are much harder to pin down. 

— li u p > u max (see equation |l6|) , the line y = u is above 
the curve |l4] at the plasma frequency. Since the curve 
|l4] goes to infinity near u z , there is an intersection near 
u z , for a frequency u > u p . The effect of cos(0) compli- 
cates this condition, which is not just u p > sub, since 
for large cos(9) even a high plasma frequency may not 
be larger than sub / \/l — cos 2 (9) cos 2 (a). 

— If in addition to the above also u x — ub < sub then 
there is also an intersection for v < u p . 

— If u p < u max , the line y — u is below the line of curve 
|l4| at u p . If, in addition, u p > sub, then there is no 
solution. This condition tells us that there are no ZM 
solutions for large cos(6>). 

— If Up < u max , but Up < sub than there is no inter- 
section in the range v p < u < u max . However, if the 
refraction index changes slowly in the vicinity of sub, 
and u z is relatively far from sub, then it is possible 
that the line y = u can intersect the curve |l4] from 
below. For such a case there is also a second solution, 
closer to u z , when the curve |lj goes to infinity. The 
condition that u z » sub is that cos(0) is small, and 
that u p is not too small. For example, if u p — 0.25^s 
then v z < 1.0308^s, and the requirement is that this 
frequency be "far" from ub / \J 1 — cos 2 (9) cos 2 (a). For 
our standard cos(a) = 0.81 this means cos(9) << 0.3. 
It is reasonable that for even smaller plasma frequen- 
cies the condition can not be met, and there is no ZM 
solution. 

— Finally, if u p < sub then there is certainly no intersec- 
tion for u < u p . 

In summary the condition for negative absorption in the 
Z-mode is either u p > sub, or cos(9) << 1 and u p /ub 
not too small. Empirically "not too small" translates as 
u p /ub > 0.25. 

An interesting phenomena for the Z-mode negative ab- 
sorption is that for 1.4 < u p /ub < 1.8 the frequency range 
of the Z-mode is between ub and 2ub- However, for this 
range of plasma frequencies the energy of the turnover [l^ 
is relatively high for s = 1, and it turns out that most 
of the curve described by ^ passes through regions of few 
electrons. The result is that the negative absorption coef- 
ficient are very small in absolute value, or that there are 
no negative absorption coefficients at all for these u p /ub 
ratios. 



6 



A. Stupp: A New Approximation Of ECM Frequencies 



6. Comparison of Approximations 

We performed numerical calculations of the absorption co- 
efficient using the power law distribution |t] with index 
S = 3, and the loss cone distribution || with cos(a) = 
0.81, cos(a — 5a) = 0.83. The calculation were performed 
with a standard magnetic field B = 360 gauss, and an am- 
bient number density which was changed to give different 
ratios of v p /vb- For every cosine of emission angle cos(0) 
to the magnetic field the largest in absolute magnitude 
negative absorption coefficient was found by scanning in 
frequency from sub to (s + \)vb- The frequency of this 
largest in absolute magnitude negative absorption coeffi- 
cient is compared with the approximation Jill and with our 
new approximation |l3| in the figures ^| to |^ 

In figure |^ the comparison is made for frequencies be- 
tween vb and 2i^b, and the figure shows that for small 
cos(0) the Dulk-Melrose approximation [IT] is very simi- 
lar to our new approximation. However, for larger cos(0), 
the new approximation is much better. The new approx- 
imation is identical with the result of the full numerical 
computation for most of the cos(#) range where there is 
negative absorption. 

In figure [?] the comparison is made for frequencies be- 
tween 2vb and 3i>b- Here the relevant cos(#) range is 
larger, and again for small cosines the Dulk-Melrose ap- 
proximation [ll] is similar to the numerical results and to 
our new approximation. However, for angles smaller than 
about 60 degrees, the Dulk-Melrose approximation begins 
to diverge, while the new approximation remains virtually 
identical with the numerical results. 

For the ratio v p = \Avb our conclusions in section ^| 
lead us to expect that negative absorption exists only for 
frequencies v > Ivb, and the numerical computations bear 
this out. Figure || shows again that for large cos(6>) the 
approximation of equation |ll] begins to diverge from the 
numerical results, while our new approximation remains 
identical to it. 

We do not show results for the Z-mode, since we are 
interested in presenting estimates for the frequencies of oc- 
currence of observable emission. The Z-mode is important 
in quenching the maser, but it is necessary to compute 
the absorption coefficients for all the modes and compare 
them to determine whether it does. 



The parameters entering the approximation are the 
angle of emission to the magnetic field 0, the loss-cone 
opening angle a, the ratio v p /vbi and the ratio v/vb- 

The new approximation can be used to define a range 
of possible frequencies of millisecond spike emission, given 
the ratio v p /vb- Or, when spike emission is detected, the 
approximation can be used to give the range of physical 
parameters in the emission region. 



7. Discussion 

We develop a new approximation for the frequencies at 
which the absorption coefficient is negative, and therefore 
the Electron Cyclotron Maser mechanism operates. This 
new approximation is given by equation |l3|, and is easy 
to compute. Our approximation gives results which are 
much more accurate than previous approximations, and 
are practically the same as the results of a full numerical 
calculation. The frequencies derived with the approxima- 
tion are within 0.01 — 0.02 vb of the numerically calculated 
frequencies. 
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Fig. 1. A loss-cone on the plane (vu,v±) is a straight line 
with angle a to the x-axis. The resonance condition de- 
fines, with the small 7 approximation, a resonance circle 
on the plane. The special circle which is tangent to the 
line of the loss-cone is the circle for the frequency where 
the largest in absolute value negative absorption occurs. 
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Fig. 2. On the plane (7,cos(0)) the solution of the reso- 
nance condition for the pitch angle (equation ||) can ap- 
pear in one of two forms. The first curve goes to +00 when 
7 approaches 1 from above, and the second curve goes to 
— 00 when 7 approaches 1. For some frequency the first 
form is tangent to the cos(</>) = cos(a) outer edge of the 
loss-cone, and this curve has the largest in absolute value 
negative absorption coefficient. 



Fig. 3. A depiction of the plot for f{y) (equation |l4|) for 
the OM, and the case vb > v p . The curve begins above 
the line y — v, and there is only one possible intersection, 



and therefore only one possible solution of equation 13 




Fig. 4. A depiction of the plot for f(y) (equation |l4|) for 
the OM, and the case vb < v v . The curve begins below 
the line y — is, and therefore there may be no solution, or 



there are two solutions of equation 13 
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Fig. 5. A depiction of the plot fo r f{v) (equation |Ti| ) for 
the ZM, and the case v p > svg / yl — cos 2 (0) cos 2 (a) . The 
curve begins below the line y = v, and therefore there may 
be one solution for v > v v or two solutions, one of them 
for v < v p . 



0.06 - 

0.05 - 
0.04 
0.03 

o.oa - 

0.01 - 



II- 
II- 



n ew 
] II 



0.78 o.no 



0.04 0.86 
cos(0) 



0.88 0.90 



0M— new app. 
OM-DM app. 
XO— new app. 
XO-DM app. 



0.4 0.5 0.6 

cos(i9) 



Fig. 7. The difference between the frequency of largest in 
absolute magnitude negative absorption coefficient found 
by numerical calculations and the frequency found by the 
Dulk-Melrose approximation [ll], and the difference be- 
tween the calculation and the new approximation [D^, as 
function of cos(0), the cosine of emission angle to the mag- 
netic field. The calculation is for the standard parame- 
ters power law index 6=3, loss cone angles cos(a) = 
0.81,cos(a — 6a) — 0.83, magnetic field of B = 360 gauss, 
and plasma frequency v v = vb- The frequencies are for 
the second harmonic between 2vb and 3isb- 



Fig. 6. The difference between the frequency of largest in 
absolute magnitude negative absorption coefficient found 
by numerical calculations and the frequency found by the 
Dulk-Melrose approximation [fl| and the difference be- 
tween the calculation and the new approximation [l3|, as 
function of cos(0), the cosine of emission angle to the mag- 
netic field. The calculation is for the standard parame- 
ters power law index 6 = 3, loss cone angles cos(a) = 
0.81, cos(a— 6a) = 0.83, magnetic field of B = 360 gauss, 
and plasma frequency v p = vb- The frequencies are for 
the first harmonic between vb and 2vb- 
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Fig. 8. The difference between the frequency of largest in 
absolute magnitude negative absorption coefficient found 
by numerical calculations and the frequency found by the 
Dulk-Melrose approximation [llj and the difference be- 
tween the calculation and the new approximation [l3|, as 
function of cos(0), the cosine of emission angle to the mag- 
netic field. The calculation is for the standard parame- 
ters power law index 6 = 3, loss cone angles cos(a) = 
0.81,cos(a— 5a) = 0.83, magnetic field of B = 360 gauss, 
and for a large plasma frequency v v = \Avb- The fre- 
quencies are for the second harmonic between 2vb and 
3v B - 



